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Abstract 

Gabor filters play an important role in many application areas for the enhancement of various 
types of images and the extraction of Gabor features. For the purpose of enhancing curved 
structures in noisy images, we introduce curved Gabor filters which locally adapt their shape 
to the direction of flow. These curved Gabor filters enable the choice of filter parameters 
which increase the smoothing power without creating artifacts in the enhanced image. In this 
paper, curved Gabor filters are applied to the curved ridge and valley structure of low-quality 
fingerprint images. First, we combine two orientation field estimation methods in order to 
obtain a more robust estimation for very noisy images. Next, curved regions are constructed 
by following the respective local orientation and they are used for estimating the local ridge 
frequency. Lastly, curved Gabor filters are defined based on curved regions and they are 
applied for the enhancement of low-quality fingerprint images. Experimental results on the 
FVC2004 databases show improvements of this approach in comparison to state-of-the-art 
enhancement methods. 

Index Terms 

Curved Gabor filter, ridge frequency estimation, curved regions, curvature, FVC2004, image 
enhancement, orientation field estimation, fingerprint recognition, verification tests, biomet- 
rics. 

1 Introduction and Related Work 

Gabor functions |18| . in the form of Gabor filters (GFs) [19] and Gabor wavelets [33j, are 
applied for a multitude of purposes in many areas of image processing and pattern recognition. 
Basically, the intentions for using GF and log-GF [16] can be grouped into two categories: 
first, GF aim at enhancing images [37] and the second common goal is to extract Gabor 
features obtained from responses of filterbanks. Typical fields of application include: 
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Texture Texture segmentation [15] and classification (35], with applications such as e.g. 
recognizing species of tropical wood [57] or classifying developmental stages of fruit flies [6T] . 

Medical and biological applications In medical imaging, GFs are applied for the en- 
hancement of structures like e.g. finger veins [60] and muscle fibers in ultrasound images [62], 
for the detection of blood vessels in retinal images [53], as well as for many other tasks like 
e.g. analyzing event-related brain activity [39], assessing osteoporosis in radiographs [42] and 
for modeling the behavior of simple cells in the mammalian visual cortex |llj . 

Optical character recognition GFs are utilized for text segmentation [27], character 
recognition |45| . font recognition |24| . and license plate recognition [5]. 

Object recognition Objects can be detected by GFs [31] . e.g. cars [36]. Moreover, GFs 
can be used for performing content-based image retrieval [3]. 

Biometrics Gabor functions play an important role in biometric recognition. They are 
employed for many physical or behavioral traits including iris |13] . face [13], facial expression 
[35] , speaker [33], speech [55], emotion recognition in speech [25] , gait [52J, handwriting [63J, 
palmprint [58], and fingerprint recognition. 

Fingerprint recognition Gabor filterbanks are used for the segmentation [JJ and quality 
estimation [31] of fingerprint images, for core point estimation [33], classification [28J and 
fingerprint matching based on Gabor features [29, 34J. GFs are also employed for generating 
synthetic fingerprints [6J. The use of GF for fingerprint image enhancement was introduced 
in [26]. 

All aforementioned applications have in common that they use straight Gabor filters, i.e. 
the x- and y-axis of the window underlying the GF are straight lines which are orthogonal. 
Having the natural curvature inherent to fingerprints in mind, we propose curved Gabor 
filters (see Figure [T]) for image enhancement |22j . A GF can be regarded as an anisotropic 
diffusion filter which smooths along the orientation and performs inverse diffusion orthogonal 
to the orientation [50]. The basic idea is to adopt the GF to the curved structure and 
smooth along the bent ridges and valleys. While this paper focuses on fingerprint image 
enhancement, curved Gabor filters might also be useful in other fields of application, e.g. for 
the enhancement of curved structures like muscle fibers, cell filaments or annual rings in tree 
discs. 

1.1 Fingerprint Image Enhancement 

Image quality [2 J has a big impact on the performance of a fingerprint recognition system (see 
e.g. [5TJ and [H]). The goal of image enhancement is to improve the overall performance by 
optimally preparing input images for later processing stages. Most systems extract minutiae 
from fingerprints [31], and the presence of noise can interfere with the extraction. As a result, 
true minutiae may be missed and false minutiae may be detected, both having a negative 
effect on the recognition rate. In order to avoid these two types of errors, image enhancement 
aims at improving the clarity of the ridge and valley structure. With special consideration 
to the typical types of noise occurring in fingerprints, an image enhancement method should 
have three important properties: 
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Figure 1: Visualization of a curved Gabor filter for a pixel near a core point. Filter parameters 
are: a x = 6.0, a y = 8.0. The underlying curved region consists of 33 parallel curved lines 
with 65 points for each line. 

• reconnect broken ridges, e.g. caused by dryness of the finger or scars; 

• separate falsely conglutinated ridges, e.g. caused by wetness of the finger or smudges; 

• preserve ridge endings and bifurcations. 

Enhancement of low quality images (occurring e.g. in all databases of FVC2004 |40j ) and 
very low quality prints like latents (e.g. NIST SD27 [20]) is still a challenge. Techniques 
based on contextual filtering are widely used for fingerprint image enhancement [JT] and a 
major difficulty lies in an automatic and reliable estimation of the local context, i.e. the local 
orientation and ridge frequency as input of the GF. Failure to correctly estimate the local 
context can lead to the creation of artifacts in the enhanced image [32J which consequently 
tends to increase the number of identification or verification errors. 

For low quality images, there is a substantial risk that an image enhancement step may 
impair the recognition performance as shown in |17j (results are cited in Table [T] of Section 
[5]). The situation is even worse for very low quality images, and current approaches focus 
on minimizing the efforts required by a human expert for manually marking information in 
images of latent prints (see [7] and [56] ) . 

The present work addresses these challenges as follows: in the next section, two state- 
of-the-art methods for orientation field estimation are combined for obtaining an estimation 
which is more robust than each individual one. In Section [3j curved regions are introduced 
and employed for achieving a reliable ridge frequency estimation. Based on the curved re- 
gions, in Section [4] curved Gabor filters are defined. In Section [5j all previously described 
methods are combined for the enhancement of low quality images from FVC2004 and per- 
formance improvements in comparison to existing methods are shown. The paper concludes 
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with a discussion of the advantages and drawbacks of this approach, as well as possible future 
directions in Section [6l 

2 Orientation Field Estimation 

In order to obtain a robust orientation field (OF) estimation for low quality images, two 
estimation methods are combined: the line sensor method [23J and the gradients based method 
[1] (with a smoothing window size of 33 pixels). The OFs are compared at each pixel. If the 
angle between both estimations is smaller than a threshold (here t = 15 ), the orientation of 
the combined OF is set to the average of the two. Otherwise, the pixel is marked as missing. 
Afterwards, all inner gaps are reconstructed and up to a radius of 16 pixels, the orientation 
of the outer proximity is extrapolated, both as described in [23]. 

Results of verification tests on all 12 databases of FVC2000 to 2004 [381 133 SO] showed 
a better performance of the combined OF applied for contextual image enhancement than 
each individual OF estimation |22j . The OF being the only parameter that was changed, 
lower equal error rates can be interpreted as an indicator that the combined OF contains 
fewer estimation errors than each of the individual estimations. Simultaneously, we regard 
the combined OF as a segmentation of the fingerprint image into foreground (endowed with 
an OF estimation) and background. 

The information fusion strategy for obtaining the combined OF was inspired by [44J. 
The two OF estimation methods can be regarded as judges or experts and the orientation 
estimation for a certain pixel as a judgment. If the angle between both estimations is greater 
than a threshold t, the judgments are considered as incoherent, and consequently not averaged. 
If an estimation method provides no estimation for a pixel, it is regarded as abstaining. 
Orientation estimations for pixels with incoherent or abstaining judges are reconstructed or 
extrapolated from pixels with coherent judgments. 

3 Ridge Frequency Estimation Using Curved Regions 

In |26j . a ridge frequency (RF) estimation method was proposed which divides a fingerprint 
image into blocks of 16 x 16 pixels, and for each block, it obtains an estimation from an 
oriented window of 32 x 16 pixels by a method called 'x-signature' which detects peaks in the 
gray-level profile. Failures to estimate a RF, e.g. caused due to presence of noise, curvature 
or minutiae, are handled by interpolation and outliers are removed by low-pass filtering. In 
our experience, this method works well for good and medium quality prints, but it encounters 
serious difficulties obtaining a useful estimation when dealing with low quality prints. 

In this section, we propose a RF estimation method following the same basic idea - to 
obtain an estimation from the gray-level profile - but which bears several improvements in 
comparison to [26j: (i) the profile is derived from a curved region which is different in shape 
and size from the oriented window of the x-signature method, (ii) we introduce an information 
criterion (IC) for the reliability of an estimation and (iii) depending on the IC, the gray-level 
profile is smoothed with a Gaussian kernel, (iv) both, minima and maxima are taken into 
account and (v) the inverse median is applied for the RF estimate. 

If the clarity of the ridge and valley structure is disturbed by noise, e.g. caused by dryness 
or wetness of the finger, an oriented window of 32 x 16 pixels may not contain a sufficient 
amount of information for a RF estimation (e.g. see Figure [3l left image) . In regions where the 
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Figure 2: A curved region (left) with 33 lines (central line in red, parallel lines start at the 
blue squares) and 65 points per line. The gray-level profile is obtained by averaging the 
interpolated gray values along each line, followed by the detection of local minima (red) and 
maxima (blue). 

ridges run almost parallel, this may be compensated by averaging over larger distances along 
the lines. However, if the ridges are curved, the enlargement of the rectangular window does 
not improve the consistency of the gray-profile, because the straight lines cut neighboring 
ridges and valleys. In order to overcome this limitation, we propose curved regions which 
adapt their shape to the local orientation. It is important to take the curvature of ridges and 
valleys into account, because about 94 % of all fingerprints belong to the classes right loop, 
whorl, left loop and tented arch [30J, so that they contain core points and therefore regions 
of high curvature. 

3.1 Curved Regions 

Let (x c ,y c ) be the center of a curved region which consists of 2p + 1 parallel curves and 
2q + 1 points along each curve. The midpoints (depicted as blue squares in Figure [2]) of the 
parallel curves are initialised by following both directions orthogonal to the orientation for p 
steps of one pixel unit, starting from the central pixel (x c ,y c ) (red square). At each step, the 
direction is adjusted, so that it is orthogonal to the local orientation. If the change between 
two consecutive local orientations is greater than a threshold, the presence of a core point 
is assumed, and the iteration is stopped. Since all x- and y-coordinates are decimal values, 
the local orientation is interpolated. Nearest neighbour and bilinear interpolation using the 
orientation of the four neighboring pixels are examined in Section [5] Starting from each of 
the 2p+l midpoints, curves are obtained by following the respective local orientation and its 
opposite direction (local orientation 8 + tt) for q steps of one pixel unit, respectively. 
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Figure 3: Comparing an oriented window (red rectangle) of 32x16 pixels as used by the 
x-signature method (left) and a curved region consisting of 33 parallel lines and 65 points 
per line. Noise can cause the x-signature method to fail, because the oriented window may 
contain an insufficient amount of information. Magnifying the window (blue) along the local 
orientation does not remedy these deficiencies in regions of curvature and would lead to an 
erroneous gray level profile. The RF estimation based on curved regions (right) overcomes 
these limitations by considering the change of local orientation, i.e. curvature. 

3.2 Curvature estimation 

As a by-product of constructing curved regions, a pixel-wise estimate of the local curvature 
is obtained using the central curve of each region (cf. the red curves in Figures H and [§. The 
estimate is computed by adding up the absolute values of differences in orientation between 
the central point of the curve and the two end points. The outcome is an estimate of the 
curvature, i.e. integrated change in orientation along a curve (here: of 65 pixel steps). For an 
illustration, see Figure |4j The curvature estimate can be useful for singular point detection, 
fingerprint alignment or as additional information at the matching stage. 

3.3 Ridge Frequency Estimation 

Gray values at the decimal coordinates of the curve points are interpolated. In this study, 
three interpolation methods are taken into account: nearest neighbor, bilinear and bicubic 
|47j (considering 1, 4 and 16 neighboring pixels for the gray value interpolation, respectively). 
The gray-level profile is produced by averaging the interpolated gray values along each curve 
(in our experiments, the minimum number of valid points is set to 50% of the points per 
line). Next, local extrema are detected and the distances between consecutive minima and 
consecutive maxima are stored. 

The RF estimate is the reciprocal of the median of the inter-extrema distances (IEDs). 
The proportion Pmaxmin of the largest IED to the smallest IED is regarded as an information 
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Figure 4: Points of high curvature are visualized by a high intensity of red. The core and 
delta points are clearly discernible in the curvature image on the right. 

parameter for the reliability of the estimation: 

max % 
__ i&IED 
Pmaxmin — '. : 

mm i 
ieiED 

Large values of Pmaxmin are considered as an indicator for the occurrence of false extrema 
in the profile (see Figure [5]). or for the absence of true extrema. Only RF estimations 
where p m axmin is below a threshold are regarded as valid (for the tests in Section [5j we used 
thr Pmaxmin < 1.5). If Pmaxmin of the gray-level profile produced by averaging along the curves 
exceeds the threshold, then, in some cases it is still possible to obtain a feasible RF estimation 
by smoothing the profile which may remove false minima and maxima, followed by a repetition 
of the estimation steps (see Figure [5]). A Gaussian with a size of 7 and a = 1.0 was applied 
in our study, and a maximum number of three smoothing iterations was performed. In an 
additional constraint we require that at least two minima and two maxima are detected and 
the RF estimation is located within an appropriate range of valid values (between g and ) • 
As a final step, the RF image is smoothed by averaging over a window of size w = 49 pixels. 
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Figure 5: The RF estimation for the depicted example is complicated by noise and a bifur- 
cation (in the bottom area of the top image). The profile obtained by averaging along the 
lines (center) contains false minima and maxima which results in the following inter-extrema 
distances (IEDs): 9, 3, 10 (between maxima), 5, 6 (between minima). Median: 6 and propor- 
tion Pmaxmin = 10/3 = 3.3. Here, false extrema are removed by one smoothing iteration with 
a Gaussian kernel (size=7, a = 1.0). IEDs for the profile at the bottom: 11, 11, 11. Median: 

11 and Pmaxmin — 1.0. 
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4 Curved Gabor Filters 



4.1 Definition 



The Gabor filter is a two-dimensional filter formed by the combination of a cosine with a 
two-dimensional Gaussian function and it has the general form: 



9(x,y,0,f,a x ,a y ) = exp 



r, \ o 



a 



cos(27r • / • xq) 



xg = x • cos 9 + y ■ sin ( 



(1) 
(2) 



ye = —x • sin 9 + y • cos 9 (3) 

In (1), the Gabor filter is centered at the origin. 9 denotes the rotation of the filter related to 
the x-axis and / the local frequency. a x and a y signify the standard deviation of the Gaussian 
function along the x- and y-axis, respectively. 

A curved Gabor filter is computed by mapping a curved region to a two-dimensional array, 
followed by a point- wise multiplication with an unrotated GF {9 = 0). The curved region 
Cij centered in (i,j) consists of 2p + 1 parallel lines and 2q + 1 points along each line. The 
corresponding array A{j contains the interpolated gray values (see right image in Figure [2]). 
The enhanced pixel E(i,j) is obtained by: 

2p+l 2q+l 

E (hj, A i,j,f(i,j)) = ^2 ^Z A ^ k ^)-9{k-p,l-q,Q,f(i tj ),<yx,<y y ) (4) 

k=0 1=0 

Finally, differences in brightness are compensated by a locally adaptive normalization 
(using the formula from [26] who proposed a global normalization as a first step before the 
OF and RF estimation, and the Gabor filtering). In our experiments, the desired mean and 
standard deviation were set to 127.5 and 100, respectively, and neighboring pixels within a 
circle of radius r = 16 were considered. 



4.2 Parameter Choice 

In the case of image enhancement by straight GFs, [26] and other authors (e.g. [H]) use 
quadratic windows of size 11x11 pixels and choices for the standard deviation of the Gaussian 
of a x = cr y = 4.0, or very similar values. We agree with their arguments that the parameter 
selection of a x and a y involves a trade-off between an ineffective filter (for small values of 
g x and <j y ) and the risk of creating artifacts in the enhanced image (for large values of a x 
and a y ). Moreover, the same reasoning holds true for the size of the window. In analogy 
to the situation during the RF estimation (see Figure [3]), enlarging a rectangular window in 
a region with curved ridge and valley flow increases the risk for introducing noise and, as a 
consequence of this, false structures into the enhanced image. 

The main advantage of curved Gabor filters is that they enable the choice of larger curved 
regions and high values for a x and a y without creating spurious features (see Figures [6] and 
[7}. In this way, curved Gabor filters have a much greater smoothing potential in comparison 
to traditional GF. For curved GFs, the only limitation is the accuracy of the OF and RF 
estimation, and no longer the filter itself. 
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Figure 6: The enhanced image 1 (second from left) of finger 10 from FVC2004 database 3 
clearly shows the capability of the curved GF to remove scars and in the enhanced image 8 
(fourth f.l.) of finger 82 from database 1 falsely conglutinated ridges are separated by the 
curved GF. Filter parameters are a x = 8.0, a y = 8.0 and the size of the curved regions is 
65 x 65 pixels. 

The authors of [63] applied a straight GF for fingerprint enhancement and proposed to use 
a circle instead of a square as the window underlying the GF in order to reduce the number 
of artifacts in the enhanced image. Similarly, we tested an ellipse with major axis 2q + 1 and 
minor axis 2p+ 1 instead of the full curved region, i.e. in Equation |4j only those interpolated 
gray values of array A^j are considered which are located within the ellipse. In our tests, 
both variants achieved similar results on the FVC2004 databases (see Table [T]). As opposed 
to [63] , the term 'circular GF' is used in [59] and [6(1] for denoting the case a x = a y . 

5 Results 

5.1 Test Setup 

Two algorithms were employed for matching the original and the enhanced gray-scale images. 
The matcher "BOZORTH3" is based on the NIST biometric image software package (NBIS) 
[53], applying MINDTCT for minutiae extraction and BOZORTH3 for template matching. 
The matcher "VeriFinger 5.0 Grayscale" is derived from the Neurotechnology VeriFinger 5.0 
SDK. 

For the verification tests, we follow the FVC protocol in order to ensure comparability 
of the results with |17j and other researchers. 2800 genuine and 4950 impostor recognition 
attempts were conducted for each of the FVC databases. Equal error rates (EERs) were 
calculated as described in [35] ■ 

5.2 Verification tests 

Curved Gabor filters were applied for enhancing the images of FVC2004 [30]. Several choices 
for a x , a y , the size of the curved region and interpolation methods were tested. EERs for 



10 




Figure 7: The detail on the left (impression 1 of finger 90 in FVC2004 database 4) is enhanced 
by Gabor filtering using rectangular windows (center) and curved regions (right). Both filters 
resort to the same orientation field estimation and the same ridge frequency estimation based 
on curved regions. Filter parameters are also identical (jp = 16, q = 32, a x = 16, a y = 32), the 
only difference between the two is the shape of window underlying the Gabor filter. Artifacts 
are created by the straight filter which may impair the recognition performance and a true 
minutia is deleted (highlighted by a red circle). 

some combinations of filter parameters are reported in Table [T] Other choices for the size 
of the curved region and the standard deviations of the Gaussian resulted in similar EERs. 
Relating to the interpolation method, only results for nearest neighbor are listed, because 
replacing it by bilinear or bicubic interpolation did not lead to a noticeable improvement in 
our tests. 

In order to compare the enhancement performance of curved Gabor filters for low quality 
images with existing enhancement methods, matcher BOZORTH3 was applied to the en- 
hanced images of FVC2004 which enables the comparison with the traditional GF proposed 
in [26], short time Fourier transform (STFT) analysis [9] and pyramid-based image filtering 
P2] (see Table 

Furthermore, in order to isolate the influence of the OF estimation and segmentation on 
the verification performance, we tested the x-signature method |26j for RF estimation and 
straight Gabor filters in combination with our OF estimation and segmentation. EERs are 
listed in the second and sixth row of Table [T] In comparison to the results of the cited 
implementation which applied an OF estimation and segmentation as described in [26], this 
led to lower EERs on DB1 and DB2, a higher EER on DB3 and a similar performance on 
DB4. In comparison to the performance on the original images, an improvement was observed 
on the first database and a deterioration on DB3 and DB4. Visual inspection of the enhanced 
images on DB3 showed that the increase of the EER was caused largely by incorrect RF 
estimates of the x-signature method. 

Moreover, we combined minutiae templates which were extracted by MINDTCT from 
images enhanced by curved Gabor filters and from images enhanced by anisotropic diffusion 
filtering. A detailed representation of this combination can be found in |21j and results are 
listed in Table [T] To the best of our knowledge, this combination performed with the lowest 
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Matcher: MINDTCT and BOZORTH3 




Enhancement method 


DB1 


DB2 


DB3 


DB4 


Original images 


14.5 


9.5 


6.2 


7.3 


Traditional Gabor filter 
Hong, Wan, Jain [26] 


(16.9) 


14.4 


7.1 


9.8 


Short time Fourier transform (STFT) analysis 
Chikkerur, Cartwright, Govindaraju [9J 


(19.1) 


11.9 


7.6 


10.9 


Pyramid-based fitering 
Ironthaier, Kollreider, Bigun |17] 


12.0 


8.2 


5.0 


7.0 


Coherence-enhancing anisotropic diffusion filtering 
Gottschlich and Schonlieb [2JJ 


10.0 


6.4 


5.0 


6.0 


Straight Gabor filters: 




x-signature based RF estimation [26] 

-1-1 -1-1 TVTTVT T — 1 A r\ A f~\ 

11 x 11, NN, F, a x = 4.0, a y = 4.0 


12.9 


9.6 


12.8 


9.5 


Curved Gabor hlters: 




01 \y Ol ATM T? /in ^r- /in 

Zl X Zl, IMlN , ill, (T x — 4.U, <Jy — 4.U 


10.2 


6.0 


4.8 


6.0 


^ v fi^ mm "P /t /in /t /in 

OO X 00, 1N1N, r , (7 X — 41.U, Oy — 41. U 


10.5 


5.9 


4.8 


6.5 


V fi^ MM TT rr /in A (\ 

OO X DO, IN IN, Hi, (Tx — 41. U, — 41. U 


9.7 


6.3 


5.1 


6.5 


Combining curved Gabor hlters and 
coherence-enhancing anisotropic dihusion hltermg [zlj: 




TV T 1 

Max rule 


9.0 


5.0 


4.2 


5.4 


Sum rule 


9.3 


4.8 


3.6 


5.2 


Template cross matching 


8.9 


4.3 


3.4 


4.9 


Matcher: VeriFinger 5.0 Grayscale 




Original images 


8.8 


6.8 


5.3 


2.0 


Curved Gabor filters: 




21 x 21, NN, E, a x = 4.0, a y = 4.0 


6.4 


3.5 


1.7 


1.3 


33 x 65, NN, F, a x = 4.0, a y = 4.0 


6.9 


3.7 


2.0 


1.4 


33 x 65, NN, E, a x = 4.0, a y = 4.0 


6.8 


3.8 


1.8 


1.7 



Table 1: EERs in % for matchers Bozorth3 and VeriFinger on the original and enhanced 
images of FVC2004 [3D]. Parentheses indicate that only a small foreground area of the 
fingerprints was useful for recognition. Results listed in the top four rows are cited from 
|17| . Parameters of the curved Gabor filters: size of the curved region, interpolation method 
(NN = nearest neighbor), considered pixels (F = full curved region, E = elliptical), standard 
deviations of Gaussian. 



EERs on the FVC2004 databases which have been achieved so far using MINDTCT and 
BOZORTH3. 

The matcher referred to as VeriFinger 5.0 Grayscale has a built-in enhancement step which 
can not be turned off, so that the results for the original images in Table [T] are obtained on 
matching images which were also enhanced (by an undisclosed procedure of the commercial 
software). Results using this matcher were included in order to show that even in the face 
of this built-in enhancement, the proposed image smoothing by curved Gabor filters leads to 
considerable improvements in verification performance. 
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6 Conclusions 



The present work describes a method for ridge frequency estimation using curved regions and 
image enhancement by curved Gabor filters. For low quality fingerprint images, in comparison 
to existing enhancement methods improvements of the matching performance were shown. 

Besides matching accuracy, speed is an important factor for fingerprint recognition sys- 
tems. Results given in Section [5] were achieved using a proof of concept implementation 
written in Java. In a first test of a GPU based implementation on a Nvidia Tesla C2070, 
computing the RF image using curved regions of size 33 x 65 pixels took about 320 ms and 
applying curved Gabor filters of size 65 x 33 pixels took about 280 ms. The RF estimation 
can be further accelerated, if an estimate is computed only e.g. for every fourth pixel hori- 
zontally and vertically instead of a pixel-wise computation. These computing times indicate 
the practicability of the presented method for on-line verification systems. 

In our opinion, the potential for further improvements of the matching performance rests 
upon a better OF estimation. The combined method delineated in Section [2] produces fewer 
erroneous estimations than each of the individual methods, but there is still room for im- 
provement. As long as OF estimation errors occur, it is necessary to choose the size of the 
curved Gabor filters and the standard deviations of the Gaussian envelope with care in order 
to balance strong image smoothing while avoiding spurious features. Future work includes 
an exploration of a locally adaptive choice of these parameters, depending on the local image 
quality, and e.g. the local reliability of the OF estimation. In addition, it will be of interest to 
apply the curved region based RF estimation and curved Gabor filters to latent fingerprints. 
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